Dark matter in the solar system II: WIMP annihilation rates in the Sun. 
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We calculate the annihilation rate of weakly interacting massive particles (WIMPs) in the Sun 
as a function of their mass and elastic scattering cross section. One byproduct of the annihilation, 
muon neutrinos, may be observed by the next generation of neutrino telescopes. Previous estimates 
of the annihilation rate assumed that any WIMPs from the Galactic dark halo that are captured in 
the Sun by elastic scattering off solar nuclei quickly reach thermal equilibrium in the Sun. We show 
that the optical depth of the Sun to WIMPs and the gravitational forces from planets both serve to 
decrease the annihilation rate below these estimates. While we find that the sensitivity of upcoming 
km^-scale neutrino telescopes to ~ 100 GeV WIMPs is virtually unchanged from previous estimates, 
the sensitivity of these experiments to ~ 10 TeV WIMPs may be an order of magnitude less than 
the standard calculations would suggest. The new estimates of the annihilation rates should guide 
future experiment design and improve the mapping from neutrino event rates to WIMP parameter 
space. 

PACS numbers: 95.35.+d,26.65.+t,95.85.Ry,96.60.Vg 



INTRODUCTION 



There is overwhelming evidence that non-baryonic 
dark matter must exist in large quantities in the uni- 
verse, yet its composition remains unknown. An inter- 
esting possibility, especially in light of anticipated data 
from the Large Hadron Collider (LHC) and the Fermi 
Gamma-Ray Space Telescope, is that dark matter con- 
sists of at least one species of elementary particle that is 
part of an extension to the Standard Model (SM) of par- 
ticle physics. Such extensions bridge the theoretical gap 
between the electroweak symmetry breaking scale (~ 10^ 
GeV) and the Planck scale 10^^ GeV), and many nat- 
urally produce WIMPs in roughly the quantities and with 
behavior consistent with observations [1] (e.g., the ncu- 
tralino in supersymmetry [2], the Kaluza-Klein photon 
in universal extra dimension (UED) models [3-5], or the 
heavy photon in Little Higgs models [6]). 

There are three main approaches to the detection and 
characterization of WIMPs. First, WIMPs may be cre- 
ated at the LHC or in future collider experiments [7]. 
Secondly, astrophysical WIMPs may be directly detected 
by measuring the energy deposited in a target during an 
interaction with a nucleus. While dark matter has not yet 
been conclusively detected in such experiments (but see 
[8, 9]), limits on WIMP-baryon cross sections are becom- 
ing ever more stringent. Current experiments have target 
masses of ~ 10 kg, and are approaching sensitivities to 
elastic spin-independent WIMP-proton cross sections of 
a^^ ^ lO"^"* cm^ [10-12]. The experiments are starting 
to cut through swaths of minimal supersymmetric stan- 
dard model (MSSM) and Little Higgs parameter space, 
but have yet to reach the cross sections possible with the 



simplest UED models [12]. In the near future, kiloton- 
sized experiments should probe spin-independent cross 
sections down to a^^ ^ lO"'"^ cm^ [13-15]. WIMPs can 
also have spin-dependent interactions with nuclei, but 
experimental limits on those are much weaker than on 
spin-independent cross sections [16-19]. 

Thirdly, it may be possible to detect particles cre- 
ated in annihilations of WIMP pairs. Many experiments 
(most obviously Fermi, but see also [20-26]) are search- 
ing for annihilation products from the Galactic center, 
the Milky Way halo, or Local Group satellite galaxies — 
the latter sites are promising because they are predicted 
to have high dark matter densities, and the annihilation 
rate goes as the square of the density. 

The Sun may also be a source of high energy neutrinos, 
arising from WIMP annihilations at its core, since its po- 
tential well is deep enough that WIMPs may be captured 
(i.e., scattered onto orbits with speeds less than the lo- 
cal escape speed from the Sun) from the Galactic dark 
halo by elastic scattering off solar nuclei. This signal 
may be observed by the next generation of large neu- 
trino telescopes [27-29].^ For the locus of models that 
are consistent with direct detection constraints and may 
produce an observable WIMP annihilation flux of neutri- 
nos, the capture rate of WIMPs in the Sun is dominated 
by spin-dependent interactions. In particular, since most 
of the mass of the Sun is in the form of hydrogen, and 
hydrogen is the only species in the Sun with a significant 
spin-dependent capture rate, neutrino telescopes should 
have very good sensitivity to cr^^ , the spin-dependent 
WIMP-proton cross section, currently the most poorly 
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^ The Earth may also accumulate WIMPs [30]. However, we have 
shown [31] that the neutrino event rates for WIMP masses and 
cross sections not yet excluded by direct detection experiments 
fall far below the thresholds for current or planned experiments. 
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constrained of all WIMP-baryon elastic scattering cross 
sections. 

The power of neutrino telescopes to detect solar 
WIMP annihilation has been demonstrated by the Super- 
Kamiokande experiment. The upper bound on the neu- 
trino flux from this experiment was used to determine a 
conservative upper limit to (7^ ^ as a function of WIMP 
mass (using a method described in [32]) that is an order 
of magnitude better than the best limit set by direct de- 
tection, assuming that the WIMP is the supersymmetric 
neutralino [33]. The next generation of neutrino tele- 
scopes will have much larger detector areas (0.1 km^ for 
Antares [27], 1 km^ for IceCube [28] and KM3NET [29], 
versus ~ 10~^ km^ for Super-Kamiokande) , albeit with 
somewhat higher energy thresholds, and should be sensi- 
tive to much lower elastic scattering cross sections than 
Super-Kamiokande (by two orders of magnitude for Ice- 
Cube [34]). 

To estimate the event rate in neutrino telescopes from 
WIMP annihilation in the Sun for a particular model, 
one needs three basic ingredients, (i) An estimate of the 
number density n(r) of WIMPs in the Sun, and hence, 
the annihilation rate Fa oc n^. (ii) An understanding of 
how the decay products, especially neutrinos and parti- 
cles that decay to neutrinos after having traveled some 
distance in the Sun, interact and propagate through the 
Sun, interplanetary space, and the Earth, (iii) A model 
of the charged-current interactions that produce muons 
in and near the detector volume and the sensitivity of 
the telescope to those muons. Thus, the event rate in a 
neutrino telescope can be described schematically by 

^ /ncutrino\ / detector \ 
$ - Pa X X . 1 

\ physics / \response / 

Ingredient (iii) is studied as part of the neutrino tele- 
scope design process. There has recently been significant 
progress in understanding ingredient (ii) [35 38] . One 
consequence of the work on this subject is that the muon 
neutrino flux at the Earth is quite uncertain for fixed 
WIMP mass and interaction cross sections. Neutrino os- 
cillations may enhance or suppress the relevant flux in 
neutrino telescopes depending on the annihilation chan- 
nel. For a given WIMP mass and annihilation rate in the 
Sun, the event rate at the Earth can vary by a factor of 
^ 10 depending on the annihilation branching ratios. 

Traditionally, ingredient (i) has bc;en trc;ated as the 
simplest of these three sub-tasks in estimating neutrino 
event rates. In the standard picture, hereafter called the 
"instantaneous thermalization model" , WIMPs are cap- 
tured in the Sun via elastic scattering off solar nuclei, 
and then quickly thermalize into a dense core at the very 
center of the Sun. The thermalization is crucial to gener- 
ating a large annihilation rate since the annihilation rate 
varies as the square of the density of the core. This model 
of rapid thermalization of all captured heavy WIMPs is 
oversimplified for two reasons. First, the thermalization 
time depends on both the WIMP-nucleon cross section 
and mass, since the cross section governs the characteris- 



tic time between scatters and the mass governs the num- 
ber of scatters required for WIMPs to reach equilibrium 
with solar nuclei. If either the elastic scattering cross sec- 
tion is small or the WIMP mass is large, thermalization 
timescales can be long. Secondly, the Sun is surrounded 
by a complement of planets. The gravitational torques 
from these planets can change or even reduce to zero — 
the frequency with which WIMPs encounter solar nuclei 
by altering the orbital paths of the WIMPs through the 
solar system. 

The goal of this paper is to investigate the validity of 
the instantaneous thermalization model and to provide 
more accurate estimates of the WIMP annihilation rate 
in the Sun. Our principal tool will be a large set of nu- 
merical simulations of the evolution of WIMPs captured 
in the Sun, including the gravitational influence of the 
planets and rescattering by solar nucleons [31, hereafter 
Paper I]. We will demonstrate that in some cases the in- 
stantaneous thermalization model overestimates the an- 
nihilation rate by a factor of ten or more. 

In Section II, we describe the standard instantaneous 
thermalization calculation of the annihilation rate of 
WIMPs in the Sun. In Section III, we briefly describe our 
simulations, which are discussed in more detail in Paper 
I. We show how the results of the simulations modify the 
standard picture of WIMP annihilation in Section IV, 
and discuss our results in Section V. 



II. THE INSTANTANEOUS THERMALIZATION 
MODEL 

The annihilation rate of WIMPs in the Sun is given by 

Ta{t) = {av)aJ d''rn^r,t), (2) 

where n(r, t) is the time-dependent number density of 
WIMPs in the Sun, and {av)a is the velocity-averaged an- 
nihilation cross section. For the remainder of this work, 
we take n(r) to be the number density of WIMPs cap- 
tured in the Sun, since in all cases where WIMP annihila- 
tion is significant, the annihilation rate due to captured 
WIMPs is many orders of magnitude greater than the 
rate due to halo dark matter streaming through the so- 
lar system. 

If the captured WIMP number density can be sepa- 
rated as 

n{r,t) = N{t)n{v), (3) 

an assumption that we show below to be valid for the 
instantaneous thermalization model, and 7T,(r) is known 
and normalized such that / n(r)d'^r = 1, then the num- 
ber N of WIMPs in the Sun can be described by 

N = C- CaN^ - CeN. (4) 
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Here, 



dC 



^/^(r,v^)/(r,v) (5) 



A 

X |v - 



Vf<Ve 



is the capture rate of WIMPs in the Sun by elastic 

scattering of halo WIMPs with distribution function / 
off solar nuclei with atomic number A and distribution 
function /a- The interaction cross section is dcr^/dO. 
WIMPs arc only considered captured if the post-scatter 
speed of the WIMP Vf is less than the local escape speed 
Vesc{r). The coefficient in the second term of Eq. (4), 
Co, is defined by 



Ta = ^CaN\ 



(6) 



If the number density can be described by Eq. (3), Ca is 
constant and has the form 



Ca = 2{av)a I d\h^{r). 



(7) 



If the number density is not in the form (3), then Ca 
will be time-dependent and it will also be necessary to 
model the time-evolution of the density profile. The 
annihilation term in Eq. (4) is twice the annihilation 
rate because most theoretically motivated WIMPs are 
self-annihilating; each annihilation event removes two 
WIMPs from the Sun. The last term in Eq. (4), CeN, is 
an evaporation rate. This term is important for very low 
mass WIMPs {m^ < 4 GeV), since Ught WIMPs may 
gain enough energy from interactions with nucleons in 
the Sun to become unbound to the solar system [39] , but 
is negligible for the range of WIMP masses (m^ > 40 
GeV) generically expected in extensions to the SM. 

If C and Ca are time-independent and Ce is negligible, 
Eq. (4) has the solution 



Nit) = y^C/Cata.nh{t/U), 



where 



(8) 



(9) 



is the timcscalc for the number of WIMPs in the Sun to 
reach equilibrium. Eq. (8) is subject to the boundary 
condition that A''(0) = at the birth of the solar system. 
The annihilation rate goes as 



Fa = ictanh'(t/te), 



(10) 



In the instantaneous thermalization model, the 
WIMPs captured in the Sun settle quickly to an equi- 
librium (i.e., the thermalization time tt <C te)- Thus, the 
WIMP density profile is separable in the sense of Eq. (3) 
and can be described by 



ne{r,t) = nc(t)e' 



(12) 



where nc(t) is determined by N(f,), m-^ is the WIMP 
mass, $(r) is the gravitational potential of the Sun, and 
T is the characteristic temperature of WIMPs in the 
Sun (typically, the core temperature of the Sun) [40 43] . 
Since the density profile is fixed as a function of time, the 
solution Eq. (8) applies, and we can express the timescale 
te in units of the age of the solar system tQ = 4.5 Gyr, 



t_Q 
t. 



C 



ll/2 r 



1030 yr- 



3 X 10- 



1/2 



(13) 



100 GeV 



3/4 



■^^ cm'^ s ^ is the annihilation cross 



Here, {av)a = 3 x 10 
section necessary to create thermal relic WIMPs at the 
abundance observed in the universe [2]. C = 10'^" yr^^ 
is the typical capture rate in the Sun by hydrogen for 
a WIMP with « 100 GeV and WIMP-proton elas- 
tic scattering cross section ap = 10"''° cm^ assuming 
that the mass density of WIMPs in the Galactic halo is 
= 0.3 GeV cm^'^ and that the local WIMP population 
is smooth and approximately described by a Maxwellian 
velocity distribution (with one-dimensional velocity dis- 
persion (7 = vq/\/2, where vq is the speed of the Local 
Standard of Rest). If ^ tha, where niA is the mass 
of the target nucleus, the capture rate goes as 



C oc p^m^aA, 



(14) 



which has the limits 



where (ta is the elastic scattering cross section. The 
dependence of the capture rate on the WIMP mass re- 
flects not only the factor of from the WIMP number 
density n-^ = p^/m^ for fixed WIMP mass density but 
also kinematic suppression since it is difficult for halo 
WIMPs to transfer enough energy to solar nuclei to be- 
come bound to the Sun. More detailed calculations of 
the capture rate can be found, e.g., in [44]. For a WIMP 
mass of = 100 GeV, the number of WIMPs in the Sun 
will have reached equilibrium unless a-p^< 10~^^ cm^ or 
(Tp^< 10~*^ cm^. Thus, the annihilation rate can be 
computed once the WIMP mass m^, spin-dependent and 
spin-independent elastic scattering cross sections, the an- 
nihilation cross section, and the local halo WIMP phase 
space density are known. 



III. THE SIMULATIONS 



^a{t) 



\C 



^C{t/te)^ = W^Caf 



if t/te 
if t/te 



> 1 
< 1 



(11) ^''^ Paper I, we described a set of simulations to deter- 

mine the lifetimes and distribution function at the Earth 



4 



of WIMPs bound to the solar system by elastic scatter- 
ing off solar nuclei. In order to understand how those 
quantities depend on WIMP mass and elastic scattering 
cross section, wc ran four sets of simulations, each with 
a different combination of WIMP mass and elastic scat- 
tering cross section (m^^ = 60 AMU, = 10~^^ cm^; 

= 60 AMU, cr^^ = lO'^^ cm^; rn^ = 150 AMU, 
cr|^ = lO-''-'' cm2; = 500 AMU, a^^^ = IQ-^^ cm^). 
In all cases, wc set = in order to simplify the in- 
terpretation of the results, but in Paper I, wc describe 
how to extrapolate the results to regimes in which spin- 
dependent scattering dominates. 

Wc modeled the halo WIMPs as having a Maxwellian 
velocity distribution in the Galactocentric frame, setting 
the one-dimensional velocity dispersion to <t = Vq/\/2. 
We used the standard solar model described in [43] to 
model the Sun. Scatters were treated in the "cold Sun" 
approximation, in which the thermal speeds of the solar 
nuclei arc neglected. This is a reasonable assumption 
since the WIMP kinetic energy K^^ ^ Ka for any nuclear 
species A, and the halo WIMP speed 3> va- 

Each simulation followed 10^ - 10^ WIMPs from the 
initial time and location of the initial scatter by a solar 
nucleus of an unbound WIMP onto a bound orbit. We 
followed the orbits using a symplectic integrator with an 
adaptive time step [45, 46], with passages through the 
Sun, close encounters with the planets, and nearly un- 
bound orbits treated as special cases. The simulations 
allowed for the possibility of additional elastic scatter- 
ing in the Sun whenever the orbits traversed the solar 
interior. In order to more easily interpret the results of 
the simulation, we used a simplified solar system con- 
sisting only of the Sun and Jupiter. Jupiter was placed 
on a circular orbit about the Sun. This restricted three- 
body problem admits a constant of motion, the Jacobi 
constant, which is a useful check on the accuracy of the 
integration code. Typically, errors in the Jacobi constant 
were oscillatory and no more than a few parts in 10^ at 
aphelion. The details of the integration scheme are given 
in Paper I. 

The integrations were terminated if the orbits met one 
of three conditions: (i) The WIMP rescattcred onto an 
orbit that was no longer Earth-crossing. Such orbits ther- 
malize quickly in the Sun, as discussed below, and hence 
were no longer relevant to our study of the WIMP dis- 
tribution hmction at the Earth, (ii) The particle was 
ejected from the solar system, (iii) The WIMP survived 
on an Earth-crossing orbit for the entire lifetime of the 
solar system, = 4.5 Gyr. Wc use the lifetime distribu- 
tions as a function of the semi-major axis of the WIMP 
orbits to construct our argument below. 

IV. INSIGHTS FROM WIMP ORBIT 
SIMULATIONS 

The simulations offer insight into two possible mecha- 
nisms that may affect the annihilation rate in the Sun: 



the finite time required to thermalize WIMPs in the Sun, 
and changes to WIMP trajectories due to gravitational 
interactions with bodies in the solar system other than 
the Sun. The former would be relevant even if the Sun 
were an isolated body, while the latter depends on the 
presence of planets in the system. We show that both 
suppress the annihilation rate by an amount that in- 
creases with increasing WIMP mass, and the suppression 
is stronger for spin-dependent than spin-independent in- 
teractions at fixed WIMP mass. 



A. Rescattering Times 

In the standard picture of WIMP annihilation, it is as- 
sumed that WIMPs thermalize with solar nuclei on very 
short timescales. In the absence of planets, the time for 
a particle to rescatter after it becomes bound to the solar 
system is 

tr ^ Px/r, (15) 

where Py^ is the WIMP orbital period and r is the optical 
depth for one passage through the Sun. The thermaliza- 
tion time is U = Xtr, where the factor A, which we 
will estimate later in this section, depends on the WIMP 
mass. If tt is longer than the age of the solar system 
^0, the particles do not thermalize in the Sun, and the 
annihilation rate of WIMPs in the Sun essentially van- 
ishes. To find the regions in WIMP parameter space for 
which tt ^ ^0, wc consider how the thcrmalization time 
depends on both the WIMP cross section and mass. 

The thermalization time depends on the WIMP-bar- 
yon cross section through the optical depth r, which is 
proportional to the cross section a. In our simulations 
CTp-^ = 10^'''' cm^ corresponds to r ~ 10~^, the exact 
value of which depends on the trajectory of the WIMP 
through the Sun. We find that for a fixed cross section, 
Tsi ^ 100 TsD (the exact value depends weakly on mass, 
for reasons described in Paper I), so that it takes a much 
higher spin-dependent cross section to reach an equiv- 
alent optical depth as a spin-independent cross section. 
This is because (i) hydrogen is the only species in the 
Sun with an appreciable spin-dependent capture proba- 
bility, and (ii) even though there are only trace amounts 
of metals (with nucleon number A > 4) in the Sun (~ 2% 
of the total mass), the spin- independent cross section per 
nucleus goes as a^^ oc (this is only moderately sup- 
pressed due to incoherent scattering). Therefore, for a 
fixed WIMP-proton cross section, the rescattering time 
is much longer if spin-dependent interactions dominate 
in the Sun rather than spin-independent interactions. 

The WIMP thermalization time depends on 
through several different effects. As increases, the 
momentum transfer in an elastic scattering event cor- 
responds to a smaller velocity change in the WIMP. It 
thus becomes more difficult to scatter halo WIMPs onto 
bound orbits. Those that are scattered onto bound orbits 
will have larger (less negative) orbital energy and hence 
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longer orbital periods. In the case of spin-dependent scat- 
tering, the median semi-major axis for m-^ = 100 GeV is 
a w 0.05 AU w 10 Rq (where Rq is the radius of the Sun 
and 1 AU w 215 i?0), a w 0.5 AU for = 1 TcV, and 
o « 5 AU for = 10 TeV. If spin-independent scat- 
tering dominates in the Sun, the median scattered semi- 

o « 0.05 AU 
10 TeV. The 



maior axis is a ss i?p for m 





for TOy = 1 TeV, and a ■ 



100 GeV, 



0.5 AU for 



median semi-major axis is well approximated by 

niA 10 



(16) 



for both spin-independent and spin-dependent scattering. 
Typically, mA ~ 10 GeV for spin-independent scattering 
and niA ~ 1 GeV for spin-dependent scattering. Since 
Kepler's third law states that oc a^/^, we find that 

3/2 

the median WIMP orbital period goes as cx . 

Once the WIMP becomes bound to the solar system, 
the number of scatters required to thermalize a WIMP 
increases with increasing mass, hence increasing X, the 
ratio of the thermalization time to the rescattering time. 
To estimate X, we consider the typical amount of energy 
lost in each scatter. In an inertial frame moving with 
the target nucleus (in the cold Sun approximation, this 
is just the heliocentric frame), the WIMP loses an energy 



Q = 



rriA 



cost 



(17) 



per scatter, where fiA ~ mAm^/ {mA+m^) is the reduced 
mass of the WIMP-nucleus system, v is the speed of the 
WIMP relative to the nucleus, and cos 6 is the center-of- 
mass scattering angle. For a high mass WIMP, ij,a ^ tua, 
so the typical WIMP energy loss is 



Q ~ niAVg 



(18) 



where Vesc is the local escape speed from the Sun (iJesc = 
618 km at the surface of the Sun, and is ~ 2x that 
value at the center of the Sun). For now, we approximate 
scatters as occurring at points in the Sun where the es- 
cape velocity is not too different from that at the surface. 
One can express the energy loss in terms of a semi-major 
axis, using 



2aa 



Q 



rriAV^ 



Solving for Uq, we find 



GMq 



rriA 



-Rf 



■0, 



(19) 

(20) 

(21) 

(22) 



to scatter, its post-scatter semi-major axis af can be de- 
scribed by 



We find 



2^f 



Q. 



(23) 



(24) 



Oj + CJg 

and hence the change in semi-major axis per scatter is 



Aa = af — tti = — 



at 



Oi -I- a„ 



(25) 



We can define a differential equation for the shrinking 
of the semi-major axis with time, 



do__ gV^ 

dt a + Gq' 



(26) 



where we have used t = P^/t as the time between scat- 
ters, and where we have expressed a in AU and t in years. 
This equation has the solution 
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1 / 3/2 3/2\ , / 1/2 1/2 



yr, (27) 



where again, all semi-major axes should be in units of 
AU. 

To estimate the X coefficient, we find X = t/r/ P^iai) 
using the time tf it takes for a WIMP with initial semi- 
major axis tti to reach a semi-major axis af ~ 2Rq, since 
the thermalization time tt is dominated by the scatters 
that occur when the orbital period is still fairly large. We 
use af = 2Rq because we consider scatters that occur 
near Rq, and because for semi-major axes much smaller 
than this, the thermal speeds of the nuclei (which we 
ignore here) become important. We restrict X to be 
strictly greater than unity, since an integer number of 
scatters is required to thermalize a WIMP. Setting o/ a 
factor of ten higher changes X by only a few percent. We 
show X as a function of o and m^/mA in Fig. 1. From 
this figure, we see that for a WIMP orbit with Oj = 1 AU, 
it takes 0(1) scatter to thermalize if m,^/mA = 100, and 
that X oc m^/mA for higher values of the mass ratio. 

We can also use Eq. (27) to estimate the lower limit on 
the elastic scattering cross section for which the median 
thermalization time of captured WIMPs tt is less than 
the age of the solar system as a function of m^/mA- 
Using Eq. (16) to describe the median semi-major axis 
of captured WIMPs, we estimate the limits to be 



10 



cm2< 



0.3 



lim 



niA J 



3/2 



1.4^ 

niA 



where we have used the fact that v^^^ = 2GMq/Rq ^ 
GMq/Rq. If a WIMP had a semi- major axis a, prior 



-0.9}, (28) 
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FIG. 1: The number of scatters X required to bring a WIMP 
with initial semi-major axis a down to af — 2Rq as a function 
of the WIMP-nucleus mass ratio. 



and 



^SD 



lim 



100 a^^ 



lim 



(29) 



If spin-indcpcndcnt scattering dominates the capture rate 
in the Sun, the two solar species that are responsible for 
most of the capture rate are helium and oxygen, so uia ^ 
10 GeV. Thus, the thermalization time exceeds the age 
of the solar system if (7^^ < 10~^^ cm^ if TOj^, = 100 GeV, 



a^' < 10-5° cm2 if = 1 TeV, and a^^ < 10"^^ cm^ if 
— 10 TeV. If spin-dependent scattering dominates in 
the Sun, the solar species that dominates the capture rate 
is hydrogen {mA ''^ 1 GeV) , so the median thermalization 



time exceeds tp, if ^10 



100 GeV, 



cr^^ < 10-*^ cm2 if = 1 TeV, and a^^ < 10-^4 cm^ 
if = 10 TeV. 



B. The Effect of Jupiter 

We now investigate how the presence of Jupiter al- 
ters the WIMP lifetimes as a function of the initial 
semi-major axes. In our simulations, the vast majority 
(~ 99.9%) of WIMPs that were initially scattered onto 
orbits with semi-major axis a < 1.5 AU rescattered on 
timescales of order t w Px/'''- It usually only required one 
scatter to reduce the semi-major axis below the threshold 
of interest because we considered spin-independent inter- 
actions only and simulated WIMPs with < 500 GeV. 
For a given WIMP mass and scattering cross section, the 
thermalization time for these WIMPs with a < 1.5 AU 
is described by Section IV A. For completeness, we note 



that there is also a long-lived population of WIMPs with 
a < 1.5 AU on secular resonances (ones which pull the 
perihelia out of the Sun for extended periods of time) that 
contribute substantially to the bound WIMP distribution 
function at the Earth, but make up only a tiny fraction 
~ 10"^ of the total population of 0.5 AU < a < 1.5 AU 
WIMPs. These WIMPs contribute negligibly to the an- 
nihilation rate. 

There are two populations of WIMPs captured onto 
bound orbits that have their lifetimes altered by gravita- 
tional torques from Jupiter. WIMPs with 1.5 AU < a < 
2.6 AU (these WIMPs do not cross Jupiter's orbit; the 
largest possible aphelion for an orbit with a < 2.6 AU is 
2x 2.6 AU = 5.2 AU, which is Jupiter's semi-major axis), 
hereafter called the "long-lived" population, have rescat- 
tering times of order tr ~ IOOP-^/t. Through a combi- 
nation of mean-motion and secular resonances, Jupiter 
pulls the perihelia of such WIMPs out of the Sun for a 
significant fraction of their lifetimes in the solar system. 
This reduces the probability of rescattering for any given 
orbital period, hence increasing the rescattering time by 
about two orders of magnitude. 

The thermalization time of this population of long- 
lived particles, tu, is dominated by the time required for 
a WIMP to rescatter to a semi-major axis a < 1.5 AU. 
To demonstrate this, we define A"i.5 to be Eq. (27) mul- 
tiplied by T / for ai — 2 AU (a typical semi-major axis 
for long-lived WIMPs) and a/ = 1.5 AU. If m^/mA — 
10'^, Xi.5 = 1, and so the time required for a WIMP 
to drop below a — 1.5 AU is ti.s = lOOXi, ^P^/t = 
IOOP^/t. Once the WIMP reaches a ~ 1.5, according 
to Fig. 1, it takes tt ~ 7P;^(1.5 AU)/t to thermalize, 
which is much less than ii.5. 

WIMPs that initially scatter onto Jupiter-crossing or- 
bits ("Jupiter-crossing population"; a > 2.6 AU) are 
ejected from the solar system on Myr timescales unless 
the optical depth in the Sun so high that the rescattering 
timescale Px/t is less than the timescale t% for torques 
from Jupiter to pull the orbital perihelia out of the Sun. 



10- 



In our simulation with m-^ — 60 AMU, 
and fjp ^ = 0, we found that ~ 75% of all Jupiter-crossing 
WIMPs were ejected before rescattering in the Sun. The 
equivalent spin-dependent cross section, if such interac- 
tions dominate in the Sun, is (t^^ ^ 10""^^ cm^. For 
the simulations with (7^^ — lO^**^ cm^ (equivalent to 
(Tp^ ^ 10^^^ cni'^), the percentage of Jupiter-crossing or- 
bits that are ejected increased to > 98%. 

For m^/mA ^ 100, it takes on average only X ^ \ 
times the rescattering time to bring a Jupiter-crossing 
orbit down to a ^ Rq. Therefore, for such mass ra- 
tios, Jupiter-crossing WIMPs will thermalize if a^^ ^ 

,2 



10-41 cm2 or a^^ > 



10 



-39 



cm . If m^/mA is much 
higher, a Jupiter-crossing WIMP is not guaranteed to 
thermalize even if it does rescatter; it may rescatter onto 
another Jupiter-crossing orbit or become a member of the 
long-lifetime population. We define X2.Q — tfr/P^ for 
tti — 4 AU and af — 2.6 AU using Eq. (27). For a mass 
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ratio of m^/mA = 10^, ~ 5, and X1.5 « 9 for the 
WIMP to go from a = 2.6 AU to a = 1.5 AU. Therefore, 
for high mass ratios {m^/niA ^ 10^), Jupiter-crossing 
WIMPs will only thcrmalizc if the clastic scattering cross 
section is significantly higher than cr^^ "^"^ 



SD , 



10-39 cm2. 



C. Mapping the Suppression in Parameter Space 

The total suppression of the annihilation rate will de- 
pend on both the WIMP mass and the WIMP-baryon 
cross section, and on whether the scattering is spin- 
dcpcndcnt or spin-indcpcndcnt. To quantify the sup- 
pression, one must determine what fraction of captured 
WIMPs belong to each population, and solve a differen- 
tial equation for the number density of WIMPs in the 
Sun given realistic thermalization times. 

We classify the types of suppression of the annihila- 
tion rate in Fig. 2 according to the description below, 
and demonstrate how to calculate the number of WIMPs 
in the Sun (and hence, the suppression with respect to 
the instantaneous thermalization model) in each regime. 
To simplify the discussion, we assume that each WIMP 
population can be described by its median thermaliza- 
tion time. A more precise calculation of the annihilation 
rates would model the thermalization time distributions 
of each WIMP population. 

(a) XPJt < t©: In this case, ah WIMPs with a < 1.5 
AU thermalize in the Sun. We define the boundary in 
cross section for this part of parameter space by consid- 
ering either the thermalization time for the median semi- 
major axis am of all captured WIMPs if am < 1-5 AU 
or the thermalization time for a = 1.5 AU if the median 
semi-major axis lies above this value. For very small 
mass ratios (m^/mA ^ 20), am is less than af — 2Rq, 
the semi-major axis we used in Section IV A to estimate 
the thermalization time. For these values, we set X = 1 
to estimate the thermalization time. 

This case breaks up into several sub-classes, depending 
on whether the Jupiter-crossing and long-lived WIMP 
populations thermalize in the Sun. 

• (al) t% > X2.ePx/T: In this case, the timescale 
for Jupiter to pull the perihelia of Jupiter-crossing 
WIMPs out of the Sun is greater than the rescat- 
tering timescale, so Jupiter-crossing WIMPs ther- 
malize in the Sun. The thermalization time cor- 
responding to X2.6 is within a factor two for the 
range of in Fig. 2 if a more careful esti- 
mate using the median Jupiter-crossing semi-major 
axis is considered. For the range of cross sec- 
tions for which Jupiter-crossing WIMPs thermalize, 
long-lived WIMPs also thermalize on timescales of 
ta = lOOXi, ^P^/t, and so the annihilation rate in 
the Sun will be unchanged from that computed in 
Section II. 

• (a2) t%. < X2.&Pxl''' and tu < Iq: In this case, 



Jupiter-crossing WIMPs will not thermalize, mean- 
ing that the capture rate of WIMPs that do ther- 
malize is C" = = {1 — f%)C, where C is the 
capture rate calculated in Eq. (5) and f% is the 
fraction of all WIMPs on Jupiter-crossing orbits. 
The capture rate may be further reduced depending 
the thermalization time of the long-lived WIMPs. 

(a21) tu < t1: If the long-lived WIMPs thermal- 
ize quickly relative to tg , the equilibrium time due 
to the reduced capture rate, then the annihilation 
rate of WIMPs in the Sun can be calculated in the 
instantaneous thermalization model of Section II, 
replacing C with in the calculation. Therefore, 
the annihilation rate will have the same form as Eq. 
(10), replacing C with C"^ . 

(a22) tu > : One can think of the capture rate of 
WIMPs in the Sun as a step-function. 



C'{t) 



{l-fu)C, t<tu 
C^, t > tu, 



(30) 



where fu is the fraction of all WIMPs initially scat- 
tered onto bound orbits that have o > 1.5 AU, and 
the capture rate reflects only those particles that 
may have thermalized on timescales < t. Solving 
the differential equation (Eq. 4) for the number of 
WIMPs in the Sun using this time-dependent cap- 
ture rate, we find 



N = \ l — tanh 



with 



{to - tu)/tl 



tanh-1 ( ^J ^No 



(31) 



iVn = 



Ca 



tanh ( V(l - MCCatu) . (32) 



• (a3) t% < X2.6Px/t and tu > None of the par- 
ticles with a > 1.5 AU will thermalize in the Sun. 

In this case, the annihilation rate has the same form 
as that shown in Eq. (10) but with a reduced cap- 
ture rate C" = (1 - fu)C. 

(b) XP^/t > tQ-. a WIMP captured onto an orbit 
with the median semi-major axis am will not thermalize 
in the Sun. 

The types of suppression have a stronger mass depen- 
dence for the range of in Fig. 2 if spin-dependent 
scattering dominates in the Sun because the typical mass 
of the nucleus on which the WIMPs scatter is a factor of 
ten smaller than for spin-independent scattering. The 
regions are also shifted up for spin-dependent interac- 
tions because r is a factor of 100 smaller than for spin- 
independent interactions for a fixed WIMP-proton cross 
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Spin-Independent Spin-Dependent 




FIG. 2: Suppression regimes, as defined in the text, as a function of WIMP mass and (left) spin-independent and (right) 
spin-dependent elastic scattering cross section. 



section. For the swath of parameter space shown in Fig. 
2, spin-independent cross sections in case (a2) are always 
in sub-case (a22) because, for fixed r, the capture rate is 
higher for spin-independent cross sections than for spin- 
dependent cross sections. Spin-dependent captures are 
more kinematically suppressed. Therefore, t'^ is system- 
atically shorter for spin-independent interactions than for 
spin-dependent interactions. 

We quantify the suppression as Fq/F^, where F° is 
the annihilation rate in the instantaneous thermaliza- 
tion model, and F^ is the annihilation rate calculated 
using the above methods. For points in WIMP parame- 
ter space for which the gravitational torques from Jupiter 
suppress the annihilation rate, the suppression depends 
on f% and /;;. In Fig. 3, we show the those crossing frac- 
tions for the initial distribution of captured WIMPs for 
both spin-independent and spin-dependent interactions. 
If spin-independent interactions dominate in the Sun, the 
fraction of captured WIMPs on Jupiter-crossing orbits is 
never high; it reaches only f% ~ 0.1 if = 10 TeV. 
The fraction of WIMPs with a > 1.5 AU is only fa ~ 0.2 
for the same WIMP mass. Therefore, we expect the sup- 
pression of the annihilation rate to be minimal is spin- 
independent interactions dominate in the Sun. 

In Fig. 4, we show the suppression as a function 
of WIMP mass for several cross sections. If o-p^= 
lO"**^ cm^, the annihilation rate is reduced by less than 
10% unless > 10 TeV. Only Jupiter-crossing WIMPs 
fail to thermalize. For such a cross section the number 
equilibrium time is short compared to the age of the solar 
system, so Fa/F° oc {1 — f%) (Eq. 11). The situation is 
virtually unchanged if (7^^ drops by two orders of magni- 
tude; the only change is that at = 10 TeV, the annihi- 



lation rate drops below equilibrium. If cTp^— 10~*^ cm^, 
the suppression increases since the long-lived WIMPs 
have thermalization lifetimes beyond Iq for all WIMP 
masses shown in Fig. 2, and the number equilibrium 
time also exceeds ^0. In this case, the annihilation rate 
is reduced by more than 10% for > 2 TeV. 

The suppression is far more pronounced if spin- 
dependent scattering dominates in the Sun. Accord- 
ing to Fig. 3, nearly all WIMPs have a > 1.5 AU if 
= 10 TeV, and 85% of captured WIMPs are on 
Jupiter-crossing orbits. Stepping through the lines in 
Fig. 4 for spin-dependent scattering, we find that for 
(Jp^— 10~^^ cm^, the suppression is linear in the capture 

rate since tQ. The annihilation rate is reduced 
by > 10% from the instantaneous thermalization model 
for ^ 1 TeV. For (Tp^= 10"''^ cm^, the suppression 
increases for 5 TeV since the number equilibrium 

timescale for the full capture rate is of order the age of 
the solar system. Thus, the annihilation rate is reduced 
from its steady-state value. The equilibrium timescale is 
increased for high masses since the total capture rate C 
decreases with WIMP mass as C oc for high m^. 

For a^^ = 10-'^^ cm^, long-lived WIMPs will have 
lifetimes exceeding the age of the solar system if ^ 
500 GeV. Moreover, the number equilibrium time of the 
WIMPs that do thermalize tg > so that the suppres- 
sion Fa/F^ (X (1 — /n)^ (see Eq. 11). The suppression 
for - 1 TeV is F^ « 0.7F°, and F^ < 0.01 r° if 
nix ^ TeV. If the cross section is much lower than 
this, 10 TeV WIMPs will not thermahze at all. 
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FIG. 3: The ratio of the capture rate of WIMPs onto orbits above a certain energy threshold to the total capture rate due to 
(left) spin-independent and (right) spin-dependent interactions in the Sun as a function of WIMP mass. 
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FIG. 4: The ratio of the estimated annihilation rate in the Sun to the annihilation rate calculated in the instantaneous 
thermalization model. In all cases, {av)a = 3 x lO"'^^ cm^. 



V. DISCUSSION 



Throughout this work, we have considered a stripped- 
down solar system with Jupiter as the only planet. We 
now ask how the above discussion is changed by the inclu- 
sion of the other planets in our solar system. The pres- 
ence of the other planets may affect the results above 



if they either (i) change the lifetimes for each popula- 
tion of WIMP (a < 1.5 AU, 1.5 AU < a < 2.6 AU, and 
a > 2.6 AU) or (ii) if they introduce new classes of WIMP 
populations. 

Given the following arguments, it is unlikely that the 
presence of other planets will affect either the lifetimes or 
the classification of WIMP orbits enough for the annihi- 
lation rate to be much different from those calculated in 
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Section IV. First, we consider Jupiter-crossing WIMPs, 
with a > 2.6 AU. Since Jupiter is by far the most mas- 
sive planet in the solar system, it largely sets the ejec- 
tion timcscalc of WIMPs from the solar system. Thus, 
the population of a > 2.6 AU WIMPs should be largely 
unaffected by the presence of the other planets. 

WIMPs with a < 2.6 AU could be affected if cither the 
resonance structure is significantly altered (recall that 
the long-lived population owes its survival to both mean- 
motion and secular resonances) or if there is a significant 
probability that most WIMPs of a given initial semi- 
major axis will experience a strong encounter with an 
inner planet. The planets other than Jupiter are not 
likely to affect the secular resonance structure for the 
highly eccentric orbits that originate in the Sun. Simula- 
tions of near-Earth asteroid orbits show that changes to 
the resonance structure due to the all non-Jupiter plan- 
ets arc only important if the orbits are initially circular, 
at low inclination, and have semi-major axes quite near 
those of either the Earth or Venus [47]. However, WIMPs 
captured in the Sun will generically have very high ec- 
centricities if they cross the orbits of the inner planets, 
so it is likely that Jupiter still dominates the resonance 
structure for most orbits. 

To estimate the importance of close planetary encoun- 
ters, we treat interactions as local and describe changes 
in the WIMP semi-major axis using a random-walk ap- 
proximation. This treatment may not be a good descrip- 
tion for the scattering of WIMPs on resonances, since 
such WIMPs may either be protected from [47] or have 
much stronger interactions with planets than predicted 
in the diffusion approximation. However, this argument 
will supply some rough interaction timescales. We con- 
sider inner planet encounters to have a significant impact 
on WIMP lifetimes only if the resulting RMS change to 
the WIMP semi-major axis {{6a)'^)/a? ~ 1. 

If we approximate the encounters to be local, the typ- 
ical change to the WIMP speed u in an inertial frame 
moving with planet P is 



5u ■ 



bu 



(33) 



for an impact parameter b. Since tlw^ WIMP orbits are 
highly eccentric, the heliocentric velocity v of the WIMP 
is nearly orthogonal to the heliocentric velocity of the 
planet Vp (since the inner planets are on nearly circular 
orbits) unless the WIMP semi-major axis a « ap, so that 



+ V 



P' 



(34) 



where v% = GMq/up is the square of the planet's orbital 
speed. Therefore, in heliocentric coordinates, the change 
to the WIMP's speed is 



6v 



GMp 
bv 



(35) 



The change in semi-major axis in each encounter is 

^2 



da 



a 



GMq 
Mp a^ 



vdv 



(36) 
(37) 



since the WIMP's energy E = -GMQ/2a = + 
*0(r-). 

The RMS change to the semi-major axis is 



{{6ar)=N{6ar, 



(38) 



where A'^ is the number of encounters with impact param- 
eter < b. To determine N, we estimate that a WIMP has 
a probability ^ (b/ap)'^ (the solid angle subtended by a 
sphere of radius b centered on the planet, as seen from 
the center of the Sun) of having an encounter with im- 
pact parameter < b during each WIMP orbital period 
Px- Therefore, if we consider the WIMP for a time t, the 
total number of encounters in this time is 



N 



t 



_b_ 
ap 



(39) 



Here, we have neglected the Coulomb logarithm In A, 
which is of order InA ~ 10. Using this factor, and com- 
bining Eqs. (37-39), we find that 



{{Saf) 



10 



Mp 



a 

ap 



t 



(40) 



Let us consider the time for {{5a)) /a^ = 1 for each 
of the inner planets. For Mercury, w 0.4 AU and 
M^/Mq k 10-^. Thus, for a WIMP with a = 1 AU, it 
takes ~ 10^^ yr for the WIMP semi-major axis to change 
significantly. For Venus and Earth, Mp/M© ~ 10~^ 
and a ^ 1 AU. The timcscalc for significant changes 
to a WIMP orbit with an initial semi-major axis of 
a = 1 AU is of order 10^" yr. Mars has a smaller 
mass (Mq, « O.IM0) and larger semi-major axis {a^ « 
1.5 AU) than the Earth, so the timcscalc for WIMP or- 
bits to change due to encounters with Mars is much 
longer than from interactions with the Earth or Venus. 
All of these timescales are longer than the age of the solar 
system, the maximum time for which wc consider WIMP 
orbits, so it is unlikely that a large number of WIMPs 
with a < 2.6 AU have thermalization properties that are 
different from those discussed in Section IV. 

Now that we have mapped out the annihilation rate 
suppression in WIMP parameter space, it is possible to 
determine how serious the suppression may be for cur- 
rent and future neutrino telescopes. For this discussion, 
we only consider changes to limits on a^^ since (i) we 
showed in Section IV that there is comparatively little 
suppression to the annihilation rate if spin-independent 
interactions dominate in the Sun and (ii) given existing 
limits on a^' , current and planned neutrino telescopes 
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cannot detect WIMP annihilation in the Sun driven by 
spin-independent scattering [e.g., 10, 11]. 

First, we consider the existing limits on (7^^ from 
the Super-Kamiokande experiment [33]. In the analy- 
sis of their data, neutrino oscillations were neglected but 
limits on the cross section were otherwise derived us- 
ing fairly conservative assumptions for a supersymmet- 
ric WIMP. The most stringent limit on the cross sec- 
tion came for ~ 100 GeV, tr^^ < IQ-^'' cm^ [33]. 
For larger masses, the limit on scales as be- 

cause (i) for such cross sections, the annihilation rate will 



WIMPs in the Sun. The expected spin-dependent cross 
section goes as [5] 



SD^ 



-2. 



have reached equilibrium, and so oc C oc cTp 
and (ii) the number of neutrinos produced in a single 
WIMP annihilation event scales approximately as 
[48]. From Fig. 2 in the previous section, it appears 
unlikely that the limit on ct^ ^ from Super-Kamiokande 
is greatly suppressed. If were right at the flux limit 
for m-)^ = 100 GeV, most WIMPs on Jupiter-crossing or- 
bits would be ejected before rescattering; however, only 
f% 10~^ are Jupiter-crossing. From Fig. 2, we see that 
for WIMPs with higher masses, even Jupiter-crossing 
WIMPs should thermalize, so the annihilation rate of 
WIMPs in the Sun should be well-described by Section 
II. 

Next-generation neutrino telescopes are anticipated to 
have ~ 100 times the sensitivity of Super-Kamiokande 

to neutrinos from WIMP annihilation in the Sun [34]. 
If we were to naively scale the limits on from 
Super-Kamiokande to km^-scale experiments, we would 
find fff ^ < 10-41 cm^ for = 100 GeV and crf^ < 
10"'^^ cm^ for = 10 TeV. In reality, from Section 
IV C, we find that the limits on (7^^ should be weaker 
by a factor of (1 — f%)^^ (Fig. 3). Alternatively, this 
means that the next-generation neutrino experiments will 
be sensitive to a smaller range of <Tp ^ than currently pre- 
dicted by the experimental collaborations. The restric- 
tion is only significant (a change in the limit by > 10%) 
if > 1 TeV. 

It should be noted that the interpretation of an ob- 
served neutrino flux, or any attempt to map the observed 
flux to an elastic scattering cross section, will need to in- 
clude the effects of the annihilation branching fractions 
and neutrino oscillation [35-38] . These latter effects may 
complicate the estimate of the annihilation rate in the 
Sun. 

One may also consider the prospects of observing 
WIMP annihilation in the Sun for specific WIMP models: 

Supersymmetry: Even in limited scans of the MSSM, 
it is apparent that a large range of (Jp^'is allowed, up 
to <t|-°~ 10"^^ cm^ [49]. Therefore, it is possible that 
MSSM neutralino annihilation may be observed with 
next-generation neutrino telescopes (e.g., Antares, Ice- 
Cube). 

Universal Extra Dimensions: In the minimal version 
of this model (only one extra dimension), cTp^ ^ CTp^, so 
we consider only the prospects for observing neutrinos if 
spin-dependent interactions dominate the capture rate of 



1.8 X 10"^^ — cm^, 



(41) 



where A is the fractional mass difference between the 
WIMP (the Kaluza-Klein photon) and the Kaluza-Klein 
quark. In order to satisfy relic abundance requirements, 
771;^. > 500 GeV, although the exact lower limit depends 
on boundary terms in the UED Lagrangian [5]. The al- 
lowed range of m^—(Tp ^parameter space straddles the 
line between cases (a2) and (a3) from Section IV C, and 
also straddles the critical number equilibrium timescale, 
te ^ ^0- Both the high WIMP escape fraction and 
the possibility that the number of WIMPs in the Sun 
is not in equilibrium drive the annihilation rate down if 
> 1 TeV. It would take a telescope with at least 
an order of magnitude more sensitivity than IceCube 
to detect even the Kaluza-Klein photon with the best 
prospects for detection. 

Little Higgs: As in UED models, a^^ » a^^ , so again 
we consider only the prospects of finding neutrinos from 
WIMPs captured in the Sun by spin-dependent inter- 
actions [50]. In Little Higgs models with T-parity, the 
natural scale for the heavy photon, the WIMP in this 
model, is m^^ < 500 GeV, and the spin-dependent cross 
section scales with mass as 



ctI-^ ~ 5 X 10 



-47 



1 TeV 



0.1 

a" 



cm 



(42) 



where A is the fractional mass difference between the 
heavy photon and the T-odd quark [6, 12, 50-53]. For 
a fiducial case of A = 0.1, ~ 5 x 10-*^ cm^ for 
= 100 GeV and cr|^ - 8 x lO^-^^ cm^ for = 
500 GeV. In the former case, the number of WIMPs 
in the Sun will marginally be in equilibrium, and the 
suppression due to WIMP populations with a > 1.5 AU 
will be negligible. One would still require ^ 10 times 
the sensitivity of IceCube to detect such a WIMP. In the 
latter case, the number of WIMPs in the Sun is small and 
far from equilibrium, but suffers little suppression due to 
WIMPs with a > 1.5 AU. However, if A is not much 
larger than the fiducial value, XP^/t > tQ for the initial 
median semi-major axis of captured WIMPs, and so the 
WIMP annihilation rate in the Sun will quickly drop to 
almost nothing. Therefore, one would expect virtually 
no neutrino signal from the Sun for a reasonable swath 
of Little Higgs parameter space. 
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